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Abstract. In this work we study the generation of water waves by an underwater sliding 
mass. The wave dynamics are assumed to fell into the shallow water regime. However, the 
characteristic wavelength of the free surface motion is generally smaller than in geophys- 
ically generated tsunamis. Thus, dispersive effects need to be taken into account. In the 
present study the fluid layer is modeled by the Peregrine system modified appropriately 
and written in conservative variables. The landslide is assumed to be a quasi-deformable 
body of mass whose trajectory is completely determined by its barycenter motion. A 
differential equation modeling the landslide motion along a curvilinear bottom is obtained 
by projecting all the forces acting on the submerged body onto a local moving coordi- 
nate system. One of the main novelties of our approach consists in taking into account 
curvature effects of the sea bed. 
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1. Introduction 

Extreme water waves can become an important hazard in coastal areas. Many geophysi- 
cal mechanisms are related with underwater earthquakes and landslides. The former genesis 
mechanism has been intensively investigated since the Tsunami Boxing Day but also be- 
fore [OkaSS, OS03, OS04, SB06, DD07, DD09, BCD+09]. The list of references is far from 
being exhaustive. In this study we focus on the latter mechanism - the underwater land- 
slides which can cause some damage in the generation region. In general, the wavelength 
of landslide generated waves is much smaller than the length of transoceanic tsunamis. 
Consequently, the dispersive effects might be important. This consideration explains why 
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we opt for a dispersive model [DKMll], which is able to simulate the propagation and 
run-up of weakly nonlinear weakly dispersive water waves on nonuniform beaches. 

Most of the landslide models which are currently used in the literature can be conven- 
tionally divided into three big categories. The first category contains the simplest models 
where the landslide shape and its trajectory are known a priori [TBCOl, THT02, LLS03, 
DNPZIO]. Another approach consists in assuming that the landslide motion is translational 
and the sliding mass follows the trajectory of its barycenter. The governing equation of 
the center of mass is obtained by projecting all the forces, acting on the slide, onto the 
horizontal direction of motion [GW99, WIGOO, DBPD09]. Finally, the third category of 
models describe the slide-water evolution as a two-layer system, the sliding mass being 
generally described by a Savage- Hutter type model [FNBB+08]. Taking into account all 
the uncertainties which exist in the modeling of the real-world events, we choose in this 
paper to study the intermediate level (i.e. the second category) which corresponds better 
to the precision of the available data. 

The present study is organized as follows. In Section 2 we briefly describe the water 
wave model we use, while Section 3 contains a more detailed presentation of the landslide 
model. The test-case considered in our study along with numerical results are presented 
in Section 4. Finally, the main conclusions of this study are outlined in Section 5. 

2. Water wave model 

The water wave model we use in this study is based on the classical system derived 
by D.H. Peregrine [Per67]. However, the original derivation assumes that the bottom is 
stationary in time, i.e. z = —d{x). Later, the bottom dynamics has been included into 
this system derivation by T. Wu [Wu81, Wu87]. In order to simulate the wave run-up, a 
conservative form of this system has to be derived. In the static bottom case it was done 
recently [DKMll]. The conservative system we use in the present study can be obtained 
in a similar way and can be written in the form: 

Ht + Q. = 0, (2.1) 

(l + -H^ — -HHxx)Qt — -H^Qxxt — -HHxQxt + + ^-^^) = gHdx + -Hdxtt, 

(2.2) 

where H{x,t) := d{x,t) + ri{x,t) is the total water depth and ri{x,t) is the free surface 
elevation below the still water level. The horizontal mass flux is denoted by Q{x,t) := 
H{x, t)u{x, t) where u{x, t) is the depth-averaged horizontal velocity variable. The bottom 
motion enters into the momentum balance equation (2.2) through the source term ^Hdxtt- 
The mass conservation equation (2.1) keeps naturally its initial form. We underline that 
the linear dispersion relation of the modified Peregrine system (2.1), (2.2) is identical with 
that the original Peregrine model [Per67] since these models differ only in nonlinear terms. 

This modified Peregrine (m-Peregrine) system (2.1), (2.2) has several advantages. First 
of all, we note that the full water wave problem is invariant under the vertical translations 
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[B082]. The asymptotic expansion method around the mean water level breaks this sym- 
metry. The introduction of conservative variables {H, Q) allows to recover this property. 
Another advantage is that the dispersive terms in the m- Peregrine system naturally vanish 
as the total water depth H along with its first derivative tend to zero. This property is 
in complete agreement with the physical behavior of water waves which become more and 
more nonlinear to the detriment of the dispersion while approaching the shore. 

In order to solve numerically the m-Peregrine system (2.1), (2.2) we choose the use of the 
finite volume method. Moreover, the run-up technique is well understood in the framework 
of Nonlinear Shallow Water equations [DPDll], which allows us to reuse this technology 
in the dispersive setting. The advective terms are discretized using the FVCF approach 
[GKLOl] with UN02 space reconstruction [H087]. The dispersive terms are treated with 
the finite differences. For the time discretization we use the Bogacki-Shampine 3rd order 
Runge-Kutta scheme with adaptive time step control. Note that on each time step we have 
to solve a tridiagonal system of linear equations in order to determine the time derivative 
Qt- We refer to [DKMll] for more details on the numerical method. 



In this section we briefly present a model of an underwater landslide motion. This process 
has to be addressed carefully since it determines the subsequent formation of water waves. 
In this study we will assume the moving mass to be a solid quasi-deformable body with a 
prescribed shape and known physical properties that preserves its mass and volume. Under 
these assumptions it is sufficient to compute the trajectory of the barycenter x = Xcit) to 
determine the motion of the whole body. In general, only uniform slopes are considered 
in the literature in conjunction with this type of landslide models [PP96, GW99, WIGOO, 
DBPD09]. However, a novel model, taking into account the bottom geometry and curvature 
effects, has been recently proposed [KSIO]. Hereafter we will follow in great lines this study. 

The static bathymetry is prescribed by a sufficiently smooth (at least of the class C^) 
and single- valued function z = —dQ{x). The landslide shape is initially prescribed by a 
localized in space function z = Co{^)- Fot example, in this study we choose the following 
shape function: 



where parameters A is the maximum height, £ is the length of the slide and xo is the initial 
position of its barycenter. Obviously, the model description given below is valid for any 
other reasonable shape. 

Since the landslide motion is translational, its shape at time t is given by the function 
z = Ci^^t) = (o{x — Xc{t)). Recall that the landslide center is located at the point with 
abscissa x = Xc{t). Then, the impermeable bottom for the water wave problem can be 
easily determined at any time by simply superposing the static and dynamic components: 



3. Landslide model 




(3.1) 



d{x, t) 



do{x) + C{x,t). 



4 



D. DUTYKH, D. MITSOTAKIS, S. BEYSEL, AND N. SHOKINA 



To simplify the subsequent presentation, we introduce the classical arc-length parametriza- 
tion, where the parameter s = s{x) is given by the following formula: 



s = L{x)= / ^l + {d',{Ordt (3.2) 

JXQ 

The function L{x) is monotonic and can be efficiently inverted to turn back to the original 
Cartesian abscissa x = L~^{s). Within this parametrization, the landslide is initially 
located at point with the curvilinear coordinate s = 0. The local tangential direction is 
denoted by r and the normal by n. 

The landslide motion is governed by the following differential equation obtained by a 
straightforward application of Newton's second law: 

= Frit), 

where m is the mass and Fr{t) is the tangential component of the forces acting on the 
moving submerged body. In order to project the forces onto the axes of local coordinate 
system, the angle 9{x) between r and Ox can be easily determined: 

9{x) = arctan((iQ(2;)) . 

Let us denote by and pe the densities of the water and sliding material correspondingly. 
If V is the volume of the slide, then the total mass m is given by this expression: 

m := {pe + c^Pu,)V, 

where is the added mass coefficient [BatOO]. A portion of the water mass has to be 
added since it is entrained by the underwater body motion. For a cylinder, for example, 
the coefficient is equal exactly to one. The volume V can be computed as 

V = W-S = W [ Coix) dx, 

where W is the landslide width in the transverse direction. The last integral can be 
computed exactly for the particular choice (3.1) of the landslide shape to give 

V = ^eAW. 

The total projected force Fj. acting on the landslide can be conventionally represented 
as a sum of two different kind of forces denoted by Fg and Fd: 

Ft = Fg + Fd, 

where Fg is the joint action of the gravity and buoyancy, while F^ is the total contribution 
of various dissipative forces (to be specified below). 

The gravity and buoyancy forces act in opposite directions and their horizontal projection 
Fg can be easily computed: 



^git) = (pe - Pw)Wg i C{x, t) sin(^(x)) dx. 
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Now, let us specify the dissipative forces. The water resistance to the motion force Fr is 
proportional to the maximal transversal section of the moving body and to the square of 
its velocity: 

1 /(is\2 
Fr = --CdPu,AWa{t)i^—) , 

here q is the resistance coefficient of the water and a{t) := sign^^^ The coefficient a{t) 
is needed to dissipate the landshde kinetic energy independently of its direction of motion. 
The friction force Ff is proportional to the normal force exerted on the body due to the 
weight: 

Ff = -Cfa{t)N{x,t). 

The normal force N{x, t) is composed of the normal components of gravity and buoyancy 
forces but also of the centrifugal force due to the variation of the bottom slope: 

N{x,t) = {pi - pu,)gW I C{x,t)cos{9{x)) dx + peW f ({x,t)K{x)(^^ dx, 

Jr Jr \dt/ 

where n{x) is the signed curvature of the bottom which can be computed by the following 
formula: 



+ 

We note that the last term vanishes for a plane bottom since k,{x) = in this particular 
case. 

Finally, if we sum up all the contributions of described above forces, we obtain the 
following second order differential equation: 

(7 + cj^l^ = (7 - l)9{h{t) - Cfa{t)m) - a{t) (cpm + ^c^a) (^)', (3.3) 
where 7 := ^ > 1 is the ratio of densities and integrals /i,2,3(^) are defined as: 

I^{t) = / ({x,t)sm{e{x)) dx, 



hit) = / Ci^^t) cos[9{x)) dx, 
Jr 

hit) = / C{^jt)K{x) dx. 
Jr 

Note also that equation (3.3) was simplified by dividing both sides by the width value W. 
In order to obtain a well-posed initial value problem, equation (3.3) has to be completed 
by two initial conditions: 

s(0) = 0, s'(0) = 0. 

In order to solve numerically equation (3.3) we employ the same Bogacki-Shampine 3rd 
order Runge-Kutta scheme as we use to solve the Boussinesq equations (2.1), (2.2). The 
integrals /i,2,3(^) are computed using the trapezoidal rule. Once the landslide trajectory 
s = s{t) is found, we use equation (3.2) to find its motion x = x{t) in the initial Cartesian 
coordinate system. 
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Bottom profile 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

x/do 



Figure 1. Bathymetry profile for parameters given in Table 1. 

4. Numerical results 

Let us consider tlie one-dimensional computational domain X = [a, 6] = [—3, 10] com- 
posed from three regions: the left and right sloping beaches surrounding a complex gen- 
eration region. Specifically, the static bathymetry function (io(x) is given by the following 
expression: 

{(Iq + 4 tan 5 ■ X, < 0, 

do + tan 5 ■ X + p{x) , < x < m, 

do + 4m tan 5 — 3 tan S ■ x, x > m, 

where the function p{x) is defined as 

p{x) = Aie-'''^"^' + A2e-*=^(^-^2)2^ 

Basically, this function represents a perturbation of the sloping bottom by two underwater 
bumps. We made this nontrivial choice in order to illustrate better the advantages of our 
landslide model, which was designed to handle general non-flat bathymetries. The values 
of all physical and numerical parameters are given in Table 1. The bottom profile for these 
parameters is depicted on Figure 1. 

The landslide motion starts from the rest position under the action of the gravity force. 
We study its motion along with the waves of the free surface up to T = 21 s. The landslide 
barycenter trajectory along with its speed and acceleration are shown in Figure 2. As it is 
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Parameter 


Value 


Voriavity acceleration, y 


i.U 


Water depth at x = 0, do 


1 n 

i.U 


Bottom slope, tan((5) 


U.iO 


uiiQeiwatei umiip ampiiLUQe, /±i 


1 n 

i.U 


Underwater bump amplitude, A2 


i.o 


Bump steepness, ki 


1 Q 
i.y 


Bump steepness, ^2 


Q n 
d.y 


Bump center position, xi 


1 
i.z 


Bump center position, X2 


Z.D 


Boundary between bottom regions, m 


4.0 


Number of control volumes, 


1000 


Slide amplitude, A 


0.3 


Length of the slide, i 


2.0 


Initial slide position, xq 


0.0 


Added mass coefficient, 


1.0 


Water drag coefficient, q 


1.0 


Friction coefficient, Cf 


tanl° 


Ratio between water and slide densities, 7 


1.5 



Table 1. Values of various parameters used in the numerical computations. 



expected, the landslide remains trapped in the second underwater bump, where it oscillates 
before stopping completely its motion. 

One of the important parameters in shallow water flows is the Froude number defined as 
the ratio between the characteristic fluid velocity to the gravity wave speed. We computed 
also this parameter along the landslide trajectory: 

Fr(t) := l"'">l 

^Jgd{xc{t),t)) 

The result is presented in Figure 3. We can see that in our case the motion remains 
subcritical during the experiment. 

We installed two synthetic wave gauges located at x = 2 (point located between two 
underwater bumps) and x = m (endpoint of the generation region and the beginning of 
the right sloping beach, m = 4 in our simulations) in order to measure the magnitude of 
generated free surface motions. These synthetic records are presented in Figure 4. Finally, 
the wave run-up is simulated numerically on the left and right beaches. The shoreline 
motion is represented in Figure 5. We can see, that the proposed scenario provides higher 
run-up values on the opposite beach to the slope where the landslide takes place. 
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3.5^ 




-0.5" 1 ' ' ' ' ' 1 

5 10 15 20 25 30 35 

t 



Figure 2. Barycenter position (blue solid line), velocity (black dashed line) 
and acceleration (blue dash-dotted line) during the landslide motion. 

5. Conclusions 

In this study we presented a novel model of a landslide motion over general curvilin- 
ear bottoms. This model takes into account the effects of bottom curvature, generally 
neglected in the literature [PP96, GW99, WIGOO, DBPD09]. Despite the inclusion of 
some new physical effects, the considered model is computationally inexpensive and can 
be potentially used in more operational context. The computed bottom motion is strongly 
coupled with a conservative Peregrine system [DKMll] which describes the propagation of 
nonlinear weakly dispersive water waves. The run-up of landslide generated waves on both 
beaches is simulated using the m-Peregrine system. The extension of this approach to three 
dimensions represents a natural perspective for future research along with investigation of 
other possible scenarios. 
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Figure 3. Local Froude number computed along the slide motion. 
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dashed line) beaches. 
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